Dissipative spin dynamics near a quantum critical point: 
Numerical Renormalization Group and Majorana diagrammatics 



S. Florens and A. Freyn 
Institut Neel, CNRS et Universite Joseph Fourier, 
25 avenue des Martyrs, BP 166, 38042 Grenoble, France 

D. Vcnturclli 

NEST, Istituto Nanoscienze, CNR and Scuola Normale Superiore, Pisa, Raly 
Institut Neel, CNRS et Universite Joseph Fourier, 
25 avenue des Martyrs, BP 166, 38042 Grenoble, France and 
International School for Advanced Studies (SISSA), Via Bonomea 265, 1-34136 Trieste, Italy 

R. Narayanan 

Department of Physics, Indian Institute of Technology Madras, Chennai 600036, India 

We provide an extensive study of the sub-ohmic spin-boson model with power law density of 
states J(u) oc u s (with < s < 1), focusing on the equilibrium dynamics of the three possible 
spin components, from very weak dissipation to the quantum critical regime. Two complementary 
methods, the bosonic Numerical Renormalization Group (NRG) and Majorana diagrammatics, are 
used to explore the physical properties in a wide range of parameters. We show that the bosonic 
self-energy is the crucial ingredient for the description of critical fluctuations, but that many-body 
vertex corrections need to be incorporated as well in order to obtain quantitative agreement of 
the diagrammatics with the numerical simulations. Our results suggest that the out-of-equilibrium 
dynamics in dissipative models beyond the Bloch-Redfield regime should be reconsidered in the 
long-time limit. Regarding also the spin-boson Hamiltonian as a toy model of quantum criticality, 
some of the insights gained here may be relevant for field theories of electrons coupled to bosons in 
higher dimensions. 



I. INTRODUCTION 

Quantum dissipative models, introduced on a phe- 
nomenological basis in many areas of physics, have 
proved to be an invaluable tool for studying the quan- 
tum mechanical friction of an exte rnal bosonic environ- 
ment on a generic two-state systemP^l While most stud- 
ies on quantum dissipative models investigated questions 
pertaining to decoherence at short time scales, their long- 
time behavior, in particular near possible quantum crit- 
ical points, is less well understood, see Refs. [3H5] for re- 
cent reviews. The motivation to investigate such regimes 
comes from the recent interest in zero temperature (quan- 
tum) phase transitions, - where anomalous low energy 
properties emerge due to the presence of quantum crit- 
ical modes. Despite their apparent simplicity, dissipa- 
tive impurity models are deeply concerned by this phe- 
nomenon, displaying often non-trivial long-time dynam- 
ics. They can be used to describe quantum criticality not 
only at the level of a single magnetic impurity diluted in 
a correlate d syst em (such as Mott insulator*^ or mag- 
netic metals^iil), but also in bulk materials themselves 
(e.g. i n qua ntum spin glasses,^ or heavy fermion com- 
pounds^^) v ia the framework of the Dynamical Mean 
Field TheoryP Although we restrict our discussion here 
to a two-level impurity model coupled to a bosonic con- 
tinuum, it is hard not to mention the host of fascinating 
physical phenomena in fermionic models, with the Kondo 
problem^ and its variations being another major exam- 
ple for dissipation (via magnetic screening) and quantum 



critical phenomena (see Ref. [3] for a recent theoretical re- 
view and Ref. [T7] for a discussion of some experimental 
realizations in quantum dots). 

Despite their quantum nature, the phase transitions in 
dissipative quantum impurity models are generally be- 
lieved to be well described by an associated long-range 
classical one-dimensional Ising modelP^ Indeed the sub- 
ohmic spin-boson model, in which a single spin 5 = 1/2 
is subjected to a magnetic field A perpendicular to the 
quantization axis, and coupled (along the quantization 
axis) to a bosonic continuum with vanishing density 
of states J{<jj) — 2Traujl~ s uj s 8(uj)9(uj c — u), obeys such 
quantum-to-classical mapping in the range 1/2 < s < 1, 
as confirmed by recent direct simulations of the quantum 
modelP^l In this case, the quantum critical behavior fol- 
lows non-trivial classical exponents of an Ising chain with 
long-range couplings, decaying as l/(r — r') 1+s in the 
imaginary time domain (r v ariabl es), and is described 
by an interacting fixed pointP^' However, the behav- 
ior in the strongly sub-ohmic case, < s < 1/2, has 
been recently debated. According to the quantum-to- 
classical mapping, the model should fall above its upper 
critical dimension and show mean-field exponents with a 
violation of hyperscaling laws. Initial evidence for a dif- 
ferent behavior, namely non-mean-field exponents above 
the upper critical dimension for the strongly sub-ohmic 
spin-boson model, came with surprise from Numerical 
Renormalization Group (NRG) calculationsP Ref. [23] 
also provided similar puzzling results from Monte Carlo 
simulations based on a truncation of the long-range in- 
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teraction (which may or may not be in the same uni- 
versality class as the spin-boson problem). At this junc- 
ture it is to be noted that these results have appeared 
recently to be controversial,^^ and at present no fi- 
nal consensus has appeared on the definite nature of this 
continuous quantum phase transition (QPT). Still, there 
are hints that the NRG displays truncation errors in the 
long-range ordered phase*2ESE21 so that the quantum- 
to-classical mapping may hold after all. We note that 
a discontinuous first order transition, as obtained from 
variational methods, 30 is likely an artifact resulting from 
a breakdown of this approximation in the quantum crit- 
ical regime, and that the localization transition in the 
spin-boson model is currently believed to be second or- 
der. 

In the present paper we will not attempt to answer all 
these unresolved issues, yet we hope to make progress on 
the quantitative understanding of the spin-boson model 
near its quantum critical point, focusing on the so-called 
delocalized regime, which corresponds to the classically 
disordered phase. Two complem entary tools will be de- 
veloped, namely NRG calculat ion d 22 * 31 * ^ and Majorana- 
fermion diagrammatic theorypESD that will be bench- 
marked against each other, both far from and close to the 
quantum critical point. Although the Majorana modes 
appear here as a technical device to reproduce the spin 
dynamics of a standard two-level system, the methodol- 
ogy developed here may be useful too in the emerging 
field of Majorana qubitsP^ Our work will also aim at an 
exhaustive study of the physical properties of the model, 
and we will examine in detail the zero-temperature equi- 
librium spin dynamics in all three possible spin direc- 
tions, which was previously not achieved to our knowl- 
edge. In this context, it is also interesting to use the 
NRG as a testbed for diagrammatic methods in fermion- 
boson models, and we will discuss how various GW-like 
schemes compare with the numerics. Interestingly, we 
will show that the phase transition is driven by a Bose- 
condcnsation of the mode mediating the interaction be- 
tween Majorana fermions, and that ladder resummation 
within the bosonic self-energy is the crucial ingredient to 
obtain the precise location of the quantum phase tran- 
sition. This allows us to establish (at two-loop order) a 
phase boundary given by the following critical dissipa- 
tion a c = (s/2 + s 2 /4)\A/u} c \ 1 - 8 + C(s 3 ), that describes 
the NRG phase diagram quite accurately for all values of 
A/uj c < 1 as long ass< 0.8. We believe that the precise 
understanding achieved here will be important to make 
further progress in the elucidation of quantum critical 
properties of the model. It is also worth stressing that the 
non- equilibrium dynamics,^ more often considered in the 
context of spin-qubitsp^ should also be deeply affected 
at long-times in case of proximity to a quantum criti- 
cal points, see e.g. Ref. 1551 for a related study. Standard 
master equation methods at the level of Bloch-Redfield 
approximations (analogous to lowest order perturbation 
theory in our equilibrium computations) should fail in 
this regime, and ought to be reconsidered seriously, pos- 



sibly at the light o f fur ther developments of both pertur- 
bative technique d 39 * 40 * and out-of-equilibrium NRGpH 

The plan of the paper is as follows. Section [TT] will 
present the general properties of the spin-boson model, 
and its solution at zero-temperature by the NRG. An 
improved broadening method 4 ^ will be used to extract 
the dynamical spin susceptibilities (the longitudinal as 
well as the transverse ones). Section III will present the 



Majorana-fermion diagrammatic method, based on per- 
turbation theory in the dissipation strength, focusing on 
the regime of very weak dissipation far from the quantum 
critical point. Then, we will extend the diagrammatics 
up to the vicinity of the quantum phase transition in Sec- 
tion |IV| Leading logarithmic corrections in the form of 
ladder diagrams will be obtained in a two-loop calcula- 
tion of the bosonic self-energy, leading to the accurate 
analytical formula for the phase boundary quoted above. 
Such renormalization effects, taking into account quan- 
titative vertex corrections, will allow to extend the Ma- 
jorana diagrammatics non-perturbatively, and match the 
numerical data both at low (critical modes) and high en- 
ergy (dissipative features) . Further directions of research 
will be proposed as a conclusion, and several appendices 
will contain some technical details, such as a derivation 
of Shiba's relation for the sub-ohmic spin-boson model, 
and an heuristic discussion of the quantum-to-classical 
mapping via the effective bosonic action. 



II. SPIN-BOSON MODEL AND ITS SOLUTION 
BY THE NRG 

A. Hamiltonian and physical aspects 

The (sub)-ohmic spin-boson Hamiltonian involves a 
three-component half-integer spin S = |<r (we set fi = I 
in what follows), with a the three Pauli matrices, and a 
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A 



A and e are magnetic fields applied to the spin (respec- 
tively orthogonal and parallel to the quantization axis), 
and the coupling constant A controls the strength of the 
dissipation to the environment. The bosonic spectrum 
(which includes by standard definition the dissipative 
coupling A) is assumed to be a power-law up to a sharp 
high-energy cutoff lo c : 



J{uj) = 2TtaLu l ( r s uj s 6(uj)0{uj c - w) 



(2) 



thus defining the dimensionless dissipation strength a — 
A 2 /(2cj 2 ), an important parameter of the model. In all 
what follows, we will consider a vanishing magnetic field 
parallel to the bosonic bath e = 0, focusing only on the 
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interplay of dissipation and spin precession, which are 
driven respectively by the parameters a and A. Tem- 
perature will also be set to zero, so that the model (T]) 
displays two distinct kinds of ground state: i) a delo- 
calized phase, with zero magnetization along the z axis, 
(S z ) = 0, in which the ground state is non-degenerate 
and is adiabatically connected to a coherent superposi- 
tion of the two spin up/down states as the dissipation 
strength a is increased from zero (keeping a finite value 
of the magnetic field A); ii) a localized phase, with a dou- 
bly degenerate ground state, where dissipation is strong 
enough to polarize the spin, leading to a non-zero mag- 
netization (5 Z ) 7^ 0, associated to a spontaneous sym- 
metry breaking. We emphasize that in the present quan- 
tum problem, long-range correlations associated to the 
slow decay in time of the bath correlations can be suffi- 
cient to generate long-range order at zero-temperature. 
When < s < 1, a quantum phase transition therefore 
occurs between the two phases at a critical value a c of 
the dissipation strength. Having in mind a perturbative 
expansion in the dissipation, we will focus in this work 
on the delocalized phase only up to the quantum critical 
point, i.e. a < a c . We also note that in both phases, a 
finite magnetization along the applied magnetic field A 
is always present, (S x ) ^ 0. 
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Several exact results are known for the model 
see the reviews in Refs. 1 1131 For s — 1 (ohmic cas~e) 
the quantum phase transition is of Kosterlitz-Thouless 
type, and is intimately connected to the ferromagnetic- 
antifcrromagnetic transition in the Kondo model. This 
connection, that can be precisely shown with the use of 
bosonizationEEH makes it possible to demonstrate that 
the critical dissipation strength is exactly a c = 1 in the 
limit A <C w c . A particularity of the ohmic case (s = 1) 
is that the crossover from coherent to incoherent spin 
dynamics (corresponding respectively to weakly damped 
and overdamped Rabi oscillations in the transverse spin 
autocorrelation functions) occurs before the onset of the 
QPT, so that the localized phase has always an inco- 
herent character, see Fig. [T] In contrast, the sub-ohmic 
model (0 < s < 1) can undergo localization before or 
after the disappearance of the coherent oscillations, de- 
pend ing o n the precise values of the microscopic param- 
os J3814U gee Fig 

[l] In particular, in the limit of small 
magnetic field A/lo c <C 1, the quantum phase transition 
occurs at small dissipation a 1, so that low-energy dy- 
namics (related to quantum criticality) and high energy 
physics (the damping of the Rabi oscillations) are gov- 
erned by different energy scales, and both can be treated 
in principle by perturbative methods. We will focus in 
this section on a Numerical renormalization group (NRG) 
solution that encompasses both weak (a <C 1) and non- 
perturbative dissipation (a ~ 0.1) regimes, before mak- 
ing comparison with analytical methods in Sections |III| 
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FIG. 1: Sketch of the generic phase diagram of the spin-boson 
model, in the ohmic (s = 1) and sub-ohmic (0 < s < 1) 
cases respectively, as a function of the dissipation parameter 
a and the dimensionless magnetic field A/w c . The super- 
ohmic (s > 1) regime shows only the crossover from coherent 
to incoherent dynamics, without a localization transition. In 
all cases, the coherent/incoherent crossover happens at com- 
parable coupling strength, namely a ~ 0.5. 



B. Bosonic NRG algorithm 

The implementation of the NRG method follows the 
standa rd pr ocedure initially introduced for fermionic 
models EU22I anc [ recently extended to bosonic Hamilto- 
niansP^ First, the spin-boson model |l]) is rewritten in a 
continuous form: 



H 



A 



de h(e) [a\ + a e 



de e a[a e (3) 



with J(e)/n = [h(e)} 2 . The bosonic bath J(e) is then 
logarithmically discretized using the Wilson parameter 
A > 1, first on the highest energy interval near the cutoff 



c ] , and then iteratively on successive decreasing 

(for n 



with i 



= A" 



-s+l, 



and IV for the weak dissipation case. 



energy windows [w„+i,w. 
strictly positive integer). This choice of discretization 
introduces the so-called z-parameter (with < z < 1) 
that is used to average over N z different realizations of 
the Wilson chairPS (this amounts to gradual changes of 
boundary conditions, allowing to obtain better statistics 
from the numerical simulations). Taking N z different z 
values, usually uniformly distributed in the ]0, 1] interval, 
allows to obtain an order N z improvement of the spec- 
tral resolution at finite energy. For very narrow spectral 
structures in the correlation functions, this method be- 
comes too expensive and and one lays recourse to im- 
proved broadening techniques, see Ref . W2\ and Sec. II C 
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The bosonic fields are then decomposed in Fourier 
modes (with p a positive or negative integer) on each 
interval to n +i < e < u n of width d n = (1 — A _1 )A _n_z 
for n > (note that do = 1 - A -1 ): 



E 



i27rpe / d n 



(4) 



This step is of course exact, and the first NRG approx- 
imation consists in neglecting all p =/= modes, keeping 
only the operator a^ n = a/ n (the full Hamiltonian is thus 
only recovered in the A — > 1 limit P^J) This leads to the 
"star" -Hamiltonian: 



H - A ^ 

-"star — 



T—\ 

71=0 V 



• a n\ + E £« a « a n (5) 



ri=0 



with the "impurity" coupling strength (the expressions 
below apply for n > only, the n = value can be 
obtained analogously) 



7, 



-- 1 'duJ(u) = 27m 1 A ' S+1) ^ A -(n+ z ) (s+ i) 



s + 1 



and the typical energy £ ra in each Wilson shell 



(6) 



IT" s + 1 1 - A~( s+2 ) 

du)U}J(uj) = : ^— —UJr.A r ' ' 



7; 



n J UJ n -\-i 



s + 2 1 - A-( s +!) 



(7) 

One then performs an exact mapping onto the so-called 
chain Hamiltonian: 



Z Z 7T 



(8) 



^ [e n blb n + t n (b\ L+l b n + blb n+1 )} 



n=0 

where T]o = J deJ(e) and the on-site energies e n and inter- 
site hoppings t n are determined numerically by a tridi- 
agonalization of the bosonic part of the star Hamilto- 
nian p)]). 

The remaining step of the NRG follows Ref. [33 adding 
successively sites into Hchain starting with n = 0. The 
site of order n roughly describes the physics for energies 
of the order w c A~™. This iterative procedure allows to 
reach exponentially small scales in a linear effort, an im- 
portant achievement for the study of Kondo physics^ or 
impurity quantum critical behaviorP It is important to 
stress that the success of the NRG relies on the expo- 
nential decrease of the chain parameters e n and t n (re- 
lated to the analogous behavior of the star parameters 
and 7„), which makes the iterative diagonalization 
reliable and stable (see however Ref. 0S] for the exten- 
sion of NRG to other classes of models). Because the 
size of the Hilbert space needs to remain finite, trunca- 
tion constraints need to be implemented. Typical cal- 
culations are usually performed with Nb = 10 bosonic 



states on each successive bosonic site of the Wilson chain 
that is added during the renormalization procedure] 33 ! 42 ! 
and N = 160 kept NRG states (truncation). Matrix 
sizes thus do not exceed 2000 x 2000 during the whole 
computation, and we proceed to N = 30 NRG itera- 
tions with A = 2 , so that a typical low energy scale 
w min = WcA~ w ~ 10~ 9 can be reached. The result- 
ing discrete energy spectra at successive NRG iterations 
are combined using the interpolation scheme proposed in 
Ref. 06] (see however the more rigorous imp lemen tation 
given by full density-matrix NRG calculations^ 4 ^) , lead- 
ing to a set of z-dependent many-body energy levels tk,z 
labelled by quantum number k. The spin-spin correlation 
function along the z-axis (i — x, y, z) is defined from the 
equilibrium spin autocorrelation functions in real time 
and their Fourier transform to real frequency: 



1 



Ci(t) = 5 <Mt), ffi (0)]) 



Ci(w) = 



do; 
2^ 



Ci(t) cos(u)t). 



(9) 
(10) 



An alternative way to introduce the spin dynamics is 
from the imaginary time spin correlation functions and 
their Fourier transform onto Matsubara frequency: 



X (r) = (Si^StiQ)) 
X (iv) = I dre WT X (r) 



X{oj + i0 + ) = x'M+«x"M 



(11) 

(12) 
(13) 



where analytic continuation was performed in the last 
equation in order to obtain the retarded spin suscepti- 
bility. These susceptibilities are simply related to the 
spin autocorrelation functions ( 10 ) by the fluctuation- 



dissipation theorem, which reads at equilibrium and for 
zero temperature (we consider the T = limit in all what 
follows) : 



Ci(u) = SignCaOxS'M. 



(14) 



We will therefore use equivalently both nomenclatures in 
the rest of the paper. There e xists an important equa- 
tion, known as Shiba's relation)- 33 * 43 ^ which connects real 
and imaginary parts of the spin susceptibility at low fre- 
quency. This exact relation is usually derived for the 
ohmic (s = 1) spin-boson model (or equivalently for the 
anisotropic Kondo model) and reads at small frequency: 



C,(w)=27raM[xi(0)] 



(15) 



We will prove in appendix[X]that this result can be gener- 
alized to the sub-ohmic case (0 < s < 1) in the following 
form 



a( w ) = j(M)bd(o)] 



(16) 



valid in the small w limit. An useful byproduct of Shiba's 
relation is that, at a magnetic ordering of the spin (cor- 
responding to the quantum critical point), the suscepti- 
bility Xz(0) diverges, so that the low-energy power law 
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behavior C z (ui) ex J(|w|) oc |w| s obeyed in the whole delo- 
calized phase should turn into a different (and diverging) 
power law. 

A last important technical step is that all these cor- 
relation functions can be computed at zero temperature 
from the raw NRG data using Lehmann's decomposition 
rule: 



k.z 



(17) 

where e 0j z is the ground state energy. This leads to a su- 
perposition of sharp (5-peaks which need to be broadened. 
This sensitive issue is examined now. 



C. Optimized broadening method 

A single NRG calculation (for a given value of z) usu- 
ally provides a set of energy levels which come by packets 
located around each Wilson shell. The energy resolution 
at scale uj n = A~"~ z is thus usually thought to be of the 
order d n = ui n — u! n +i = (1 — A _1 )w„, and degrades at 
higher energy. In order to generate smooth NRG spectra, 
the delta-peaks in the Lehmann formula (171 are there- 
fore usually broadened^ at energy u> n on a scale of the 
same order (the broadening scale is typically buj n , with 
b ~ 0.7) using the substitution: 



S (\uj\ - u n ) 



-[i°g(MA^)7f>] 2 



(18) 



As a result, spectral features at a given energy to that 
are sharper than ui itself cannot be well resolved, and 
will come out overbroadened. Some relative degree of im- 
provement can be achieved by combining N z NRG runs 
with different values of the z parameter (z-averaging pro- 
cedure, or z-tridP^J). Although the mathematical justi- 
fication for this procedure is still unclear, in the optimal 
case the resulting energy packets turn out to be uniformly 
distributed if the set of z-values is also chosen uniformly, 
and this allows a decrease of the broadening parameter 
to b ~ 0.7/N z . In the case of very sharp features, this 
method requires a large number of NRG runs, and will 
be prohibitive. It may also become problematic in re- 
gions where the z-averaging does not provide a uniform 
distribution of states (this may occur from the accidental 
disappearance of some NRG states), in which case uncon- 
trolled oscillations of period A can be generated (situa- 
tion of underbroadening) . These problems have led au- 
thors to speculate that quantitative NRG spectral func- 
tions can be extracted only in the continuum limit B2HH] 
either with A — > 1 or N z — > +oo (only the former is 
mathematically sound, but the NRG algorithm cannot 
be managed anymore as scale separation breaks down). 
Both limitations of the z-averaging can however be lifted 
thanks to a simple remark made by two of us in a previ- 
ous publication: 42 for a single NRG run, the levels within 



a given logarithmic energy window are not uniformly dis- 
tributed, but rather tend to bunch together close to sharp 
resonances. This effect becomes quite obvious when the 
width of such resonance becomes extremely small, see 
Fig. [2] but is always present in all NRG calculations. 
Figure [2] shows indeed the raw NRG spectra for one of 
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FIG. 2: (Color online) The transverse (along the bath) 
spin susceptibility C z (oS) of the sub-ohmic spin-boson model 
at s = 0.1, A = 0.1cj c and a = 0.000125. Raw data 
|(0, z\a z \k, z)| 2 /[2iV z (e fciZ — eo,z)] are given as circles for N z = 
20 combined NRG calculations, solid line is the perturbative 
result, and the three dashed lines are the various NRG broad- 
enings discussed in the text. All frequency-dependent plots 
are given in units of u c — 1. 

the spin susceptibilities of the spin-boson model, com- 
puted with N z — 20 interleaved z-averaged NRG runs 
(dots) and smoothed using two standard NRG broad- 
enings (dash-dotted curves) with respective parameters 
b = 0.7 and b = 0.7/N z = 0.035. The largest broad- 
ening clearly fails to produce the expected peak (as a 
testbed, an accurate analytical expression is provided by 
the full line, see Sec. Ill for further details), while the 
smallest broadening signals the peak, but overestimates 
its width by more than one order of magnitude! Large 
scale z-averaging with N z ~ 10 4 would probably allow 
to resolve the correct spectral structure, but is unreason- 
ably expensive. The key observation^! is that the NRG 
eigenvalues do not always come in packets uniformly dis- 
tributed within the Wilson shell, but on the contrary tend 
to cluster close to resonances, as clearly seen on the raw 
data (dots) in Fig. [2j This provides the missing informa- 
tion which is needed in order to achieve a good spectral 
resolution with a limited numerical cost. The main strat- 
egy to be used is that the broadening parameter b must 
become energy- dependent, which we previously called the 
"b-trick" E2 However, the generic implementation of this 
simple idea is not fully clarified: we present here a simple 
scheme that seems to be quite efficient, but we empha- 
size that a more generic and robust method is still to 
be found. The required algorithm must indeed adapt 
the broadening parameter b{uj) in Eq. ( 18 ) to the fre- 



quency dependence of the density of highest weight NRG 
peaks. Our approach, which draws inspiration from the 
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NRG logarithmic discretization, is to extract b(u>) from 
the logarithmic derivative of the integrated spectrum up 
to frequency u>: 



&(«) 



bo 
2 



dlogw 



dlog/r C 
dlogo; 



(19) 

where q ~ 1 is a regularization parameter whose precise 
value does not matter much to the final result, and 60 pro- 
vides the typical broadening at low and high frequencies 
(far from the atomic resonances). Note that we have used 
here two different frequency sweeps, one from uj = and 
one from ui = +00, in order to treat on an equal basis low 
and high frequency tails (this procedure is not satisfac- 
tory in cases where several resonances are present). Be- 



cause the actual NRG data is fully discrete, see Eq. (17 1, 
we compute b(u>) recursively using Eq. (|19[) on the broad- 



ened NRG spectra. This procedure converges after few 
iterations to the results displayed on Fig. [2] It is clearly 
remarkable that such a good agreement can be reached 
with so little numerical effort, demonstrating that the 
raw NRG data encode much more information than pre- 
viously thought (see Ref. 42 for further details). One im- 
portant advantage of the b-trick is that b(u>) takes very 
small values near resonances, greatly enhancing the reso- 
lution, while keeping large values away from the peak, 
thus avoiding the usual underbroadened NRG oscilla- 
tions. The main tuning parameter in our method is the 
low-frequency broadening &o m Eq. (19), which is deter- 
mined by a simple convergence method. We finally em- 
phasize that some degree of z-averaging is always needed 
in order to generate smooth spectra with the b-trick, but 
in our experience z-averaging alone is never sufficient to 
obtain accurate results near sharp finite-frequency reso- 
nances. 



in Fig. [3] (the plot also includes a comparison to lowest or- 
der pertubation theory, introduced in Section III B ) . The 
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D. NRG results for the dynamical susceptibilities 
in the three spin directions 

We present here the NRG data for a fixed value of the 
magnetic field A/uj c — 0.1, varying the dissipation from 
values a <C a c (far from the quantum critical point) up 
to a c . We choose three different values of the bath expo- 
nent s = 0.1, 0.5, 1.0, which span the complete generality 
of the model. We remark that the case s = is special, as 
the spin is localized for all a > at zero temperature^, 
and will not be considered here, as we aim at understand- 
ing the delocalized phase of the model only. We will also 
not examine the situation of s > 1 , which shows for weak 
dissipation similar features as the s = 1 case, and does 
not present a phase transition at larger coupling. We 
stress already that depending whether a c < 1 or not, 
the transition will or won't take place at weak coupling, 
and this can lead to different physical pictures. 

Let us consider first the case s = 0.1 and the spin sus- 
ceptibilities along each spin direction (i = x, y, z) shown 



FIG. 3: Spin susceptibilities C z {uj), C v {uj) and C x {uj) (top to 
bottom panels) for A/uj c = 0.1, s — 0.1 and increasing dissi- 
pation strength a = 0.001, 0.004, 0.007 (bottom to top). The 
critical point is located at a c — 0.007. The NRG plots were 
generated with N z — 20 interleaved z-averages and used the 
b-trick discussed in Sec. |II C| while the perturbative expres- 
sions (G°W° scheme) are given in Eqs. p6p7| . Note that 
the upper panel has a different range in uj than the two other 
panels. 

transverse (orthogonal to the field and parallel t o the 
bath) spin susceptibility C z (cj), already considered, 1 41 1 42 1 
shows two interesting features. First, the resonance at 
uj = A broadened by the bath is always sharply defined, 
because here a c ~ 0.007 is so small that one remains 
in the perturbative regime for all a < a c values regard- 
ing this high energy features. This is consistent with 
the near perfect agreement with the perturbative result 
in this range of frequencies, and weakly damped oscilla- 
tion should occur in the real-time spin dynamics. The 
low-energy part of this spectral function is yet more in- 
teresting. For a <C a c , a slow decay with C z (uj) oc \ui\ s 
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is observed, also consistent with the perturbative result. 
Increasing a, a crossover at a scale uj* is seen between 
the |w| s behavior at u < w* and a different power law 
\oj\~ s at « w « A. At the quantum critical point 
a = a c , uj* vanishes and a complete |w|~ s divergence 
is realized throughout. Clearly, the perturbative expres- 
sion fails to reproduce this behavior. We stress here that 
even values relatively far from the critical point, such as 
a = 0.004, show large quantitative deviation from lowest 
order perturbation theory at low frequency, despite the 
good agreement on the coherent peak. This shows that 
the range of validity of perturbative methods (such as 
Bloch-Redfield) , although possibly accurate at short time 
scales, is very limited in the long time limit. Clearly the 
high energy dissipation mechanism does not care for the 
complex low-energy behavior of the spin dynamics, even 
at the quantum critical point PUSH Xurning to the sec- 
ond transverse susceptibility C y (u>) (orthogonal to both 
the magnetic field and the spin-bath term), we observe 
similar behavior as in C z (ui) for the resonant peak, and 
good agreement with perturbation theory. The low en- 
ergy part of the spectrum is however much less dramatic, 
because critical modes, related to fluctuations driven by 
the bath, pertain mainly to the z-component of the spin. 
What is actually going on is a very mild crossover from 
uj 2+s at w « w* to lu 2 ~ s at uj* < uj < A (the latter be- 
havior includes the whole range of low frequencies at the 
critical point), due to the small value of s = 0.1. Agree- 
ment with perturbation theory appears thus much better, 
although small deviations, due to the tiny changes of the 
power law exponent, can be seen at low frequency. These 
precise (and exact) values for the various power laws will 
be demonstrated in Section \LU\ Finally, the longitudinal 
(parallel to the bath) spin susceptibility C x (oS) presents 
quite different features. Clearly in the absence of dissi- 
pation (a — 0), a S(uj) peak at zero frequency occurs, 
in contrast to the S(oj — A) peak at the magnetic field 
frequency displayed by both C y (uj) and C z {uj). Thus no 
sharp resonance is generated at small dissipation, and 
rather a broad shoulder emerges, quite reminiscent of 
the T-matrix of the Kondo model in a magnetic fieldPSl 
This structure is well accounted for by the perturbative 
calculation for A < u, however, stark discrepancy is seen 
at low frequency for \uj\ < A: the NRG data present a 
tail of low energy modes, while perturbation theory pro- 
vides a gap in this range. This new feature given by the 
NRG will be elucidated in Section IIII CI in connection 



with multiparticle effects in the diagrammatics, that are 
neglected in the current lowest order calculation. 

We now turn to the intermediate value of s = 0.5, giv- 
ing rise to a somewhat larger critical dissipation strength 
a c ~ 0.105, see Fig. [4} In that case, the smallest chosen 
value of a — 0.02 still lies in the perturbative regime, 
so that again near perfect agreement with perturbation 
theory can be checked. However, the approach to the 
quantum critical point shows larger quantitative devia- 
tions, both at high and low energy. Clearly, the coher- 
ent peak is more strongly smeared in the NRG for the 




0.001 



FIG. 4: Spin susceptibilities C z (ui), C y (uj) and C x (uj) (top 
to bottom panels) for A/cj c = 0.1, s = 0.5 and increasing 
dissipation strength a — 0.02,0.06,0.1 (bottom to top). The 
critical point is located at a c — 0.105. The NRG plots were 
generated with N z — 10 interleaved z-averages and used the b- 
trick discussed in Sec. |II C| while the pert urbative expressions 
(G°W° scheme) are given in Eqs. |36|37j . Note that the upper 
panel has a different range in u than the two other panels. 



larger shown values a = 0.06,0.1 compared to perturba- 
tion theory, due to the now relatively important magni- 
tude of a. Nevertheless, the peak structure survives up 
to the critical point, so that moderately damped oscilla- 
tions should show up in the time domain, even close to 
a c . Regarding the low energy part of the spectrum, the 
same crossover from uj s to uj~ s behavior (for uj <C uj* and 
a;* <C uj <C A respectively) is seen in C z {lu), now more 
clearly defined due to the larger value of s = 0.5. Again 
perturbation theory fails in this regime, and this is now 
more easily seen on C v (ui) as well, which presents respec- 



tively a crossover from uj 2+s to 



in the frequency 



dependency. We note also a similar enhancement of the 
low energy tails for the last component C x {uj) below the 
shoulder at uj < A. 



We finally consider the so-called ohmic case s = 1, 
for which a c ~ 1, clearly out of the perturbative range, 
see Fig. [5] Again, small values of a = 0.1 present well- 




FIG. 5: Spin susceptibilities C z (uj), C y (uj) and C x (u>) (top to 
bottom panels) for A/w c = 0.1, s = 1 and increasing dissipa- 
tion strength a = 0.1,0.5,0.9 (bottom to top). The critical 
point is located at a c — 1. The NRG plots were generated 
with iV z = 3 interl eaved z-averages and used the b-trick dis- 
cussed in Sec. II C while the perturbative expressions (G°W° 
scheme) are given in Eqs. (36|37) . 



crossover from lu 2+s = lu 3 to ui 2 ~ s = u> power laws for 
u <C uj* and a;* <C uj <C A respectively. Finally C z {ui) 
presents also a similar crossover in the low-energy tails, 
and a broad shoulder at the magnetic field value (for 
weak dissipation). 

We now turn to the analytical description of this 
physics, and show that both qualitative and quantitative 
aspects of the numerical data can be well understood on 
the basis of the Majorana-fermion diagrammatic method. 



III. MAJORANA DIAGRAMMATICS 

A. Spin representation and correlation functions 

The idea we henceforth develop is to use a Majorana 
descriptio n 134 * 35 ^ for the impurity quantum spin S = \a, 
in order to capture quantitatively the spin fluctuations in 
the whole delocalized phase up to the quantum critical 
point using a weak dissipation expansion (understand- 
ing the quantum phase transition coming in from the 
localized phase is a more delicate issue, that we will 
not attempt here). Technically, our methodology in- 
volves the introduction of a triplet of real fermions rjj 
(with j — x,y,z) that satisfy the anticommutation rela- 
tions {rjj, r/k} = Sjk, leading to the faithful representa- 
tion Sj = —(i/2) Yli-i CjkiVkVi- We emphasize that Ma- 
jorana quoits^ involve a doublet of real fermions, while 
the mathematical decomposition of a two-level qubits re- 
quires a triplet of real fermions. The main motivation 
behind considering the Majorana representation instead 
of more standard ones (Abrikosov or Popov fermions, 
Holstein-Primakoff or Schwinger bosons) is that one can 
alternatively write S — xVj where x = SiilxVyVz is a 
fermionic operator which commutes with the Hamilto- 
nian ([I]). This implies the crucial property that the spin 
susceptibilities Eq. ( 10 ) are obtained as simple identities 



involving standard one-particle Majorana Green's func- 
tion (here at zero temperature): 



Ci(uj) = -~lmGi(ui), 

7T 



(20) 



defined resonance and reasonable agreement with per- 
turbation theory, although some deviations are visible. 
However, all values of a > 0.5 show a complete disap- 
pearance of the coherent peak at uj = A, leaving only 
the crossover scale uj*, which separates regimes where 
C z (ui) oc ui for u> <C to C z (u>) oc a;" 1 (up to loga- 
rithmic corrections) for w* C w < w c . In that case, 
oscillations are completely overdamped, and coherence is 
lost already well before the quantum critical point a c = 1 
where the spin becomes logarithmically free (this point 
corresponds precisely to the ferromagnetic scaling limit 
in the associated Kondo model). Again, C y (ui) shows a 
similar behavior as observed for smaller s values, with a 



where Gi(ui) is the retarded real- frequency Majorana- 
fermion Green's function, obtained from the Fourier 
transformation of the time-ordered Green's function 
Gi(t) = —i6(t)({rii(t) , rji(0)}) . This key relation between 
Ci and Gi is in great contrast to all other spin represen- 
tations where the spin correlation functions are given in 
terms of four-fermion correlators, making the diagram- 
matics much more cumbersome. Therefore, thanks to 



Eq. ( 20 ) , both the dissipative features at high energy and 



the low-energy quantum dynamics close to the quantum 
critical point can be simply encoded via Dyson equation 
by a Majorana-fermion self-energy. To obtain these im- 
portant relations, we need to reexpress the spin-boson 
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\fOL = 



FIG. 6: Bare vertex of the spin-boson model in the Majorana 
language. 



model Eq. ([IJ in the Majorana language: 
H = -iArj y r) z - iXr) x r) y }](a] + a^) + ^ w ' a ] a r ( 21 ) 



The "diagonal" part given by the magnetic field is read- 
ily diagonalized into "free" (or bare) Majorana Green's 
functions: 



G» 



1 



(lj + i0+) 2 - A 2 
(lj + i0+) 2 - A 2 



(22) 
(23) 
(24) 



which are readily interpreted by a spin precession in the 
y — z plane around the magnetic field pointing in the 
x direction. Interaction between the Majorana fermions 
is generated by the coupling to the bosonic bath, which 
can be expressed as a bare vertex involving the rj x and rj y 
field, see Fig. [6j This coupling results in three Majorana 
self-energies Sj(cj), but clearly the z-component vanishes 
altogether, S z (w) = 0, due to the form of the interaction. 
Inverting the Majorana propagator (written in matrix 
form): 



G = 



w-E x (w) 
lj — iA 
-iA lj 



(25) 



provides Dyson's equation for the full retarded Majorana 
Green's functions that encode the complete spin dynam- 
ics: 



G x (lj) 

Gy(Lj) 

G z (lj) 



1 



LJ - E x (u>) 
LJ 

lj[lj- - A 2 

LJ — Sy(w) 

lj[lj — T, y (u>)] — A 2 



(26) 
(27) 

(28) 



The bare Green's functions ( 22p4 ) are correctly repro- 
duced in the zero-dissipation limit, replacing £j(w) by a 
vanishing imaginary part — i0 + . 

Because the interaction term in the Hamiltonian 
Eq. ( |21[ ) involves a particular combination of the bosonic 
modes, it is useful to define the full bosonic propa- 
gator G&(t) = -i0(t) {[$(t), $(0)]) of the "local" field 



$ = —i-j= y^„-(<zj + a-). This allows us to introduce a 
bosonic self-energy E$(w) for the full Green's function 



G*(w) 



1 



G%(lj) 
G°(w)E*(w) 



in terms of the bare bosonic propagator 

A 2 A 2 



G%(lj) 



3 



a 



LJ - 



(29) 



(30) 



Using the convention chosen to normalize $, the bare 
vertex takes simply the value yfa. (see Fig. [6]), and the 
bare bosonic Green's function can be expressed through 
the bosonic spectral density J(lj) given by Eq. 



G%(lj) = - 



de J(|e|)Sign(e) 



-2lj 



l-s 



7r lj + i0 + — e 

e| s Sign(e) 



de 



LJ - 



i0+ - e 



(31) 



a quantity which is indeed independent of the dissipa- 
tion strength a, now encoded (by convention) only in 
the bare vertex. We thus find an explicit expression for 
the imaginary part of the bosonic Green's function: 



lmG%(Lj) = 27raj c Sign(w) 



0(^ c 



■ LJ 



(32) 



The real part of the integral (31 ) can be computed ana- 
lytically from Kramers-Kronig's relation only in the low- 
energy limit \lj\ <C lj c : 



4fd, 



1ZeG%(Lj) — — - — 27rw c cotan 





LJ 


(t) 


LJ C 



(33) 



although the complete formula (31 ) will be used for later 
numerical computations. The low-energy expression (33) 
breaks down for 8=1, where additional logarithmic cor- 
rections arise (hinting at the known connection to the 
Kondo modeP). 

The exact resummation of perturbation theory re- 
quires the accurate knowledge of the full fermionic and 
bosonic self-energies, as well as the full three-particle ver- 
tex r (a quantity that depends on three independent ex- 
ternal frequencies). The general diagrammatics is given 
formally in Fig. [7] We note that replacement of the bare 
propagators by the full ones without considering together 
the vertex corrections must be done in general with care, 
and can lead to serious pathologies, see Ref. [S21 An ex- 
ception to this rule is given in Sec. Ill C and the pre- 



cise role of the interaction vertex will be elucidated in 
Sec. lTVEl 



B. Lowest order perturbative self-energies in the 
Majorana diagrammatics 

For very weak dissipation a«l and far from the quan- 
tum critical point a <C a c , the spin dynamics should in 
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FIG. 7: Exact formal diagrammatic expression for the 
fermionic and bosonic self-energies in terms of the full (exact) 
fermionic and bosonic Green's functions (d and G$ respec- 
tively) and the full (exact) three-particle vertex F. 
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FIG. 8: Lowest order perturbative Majorana self-energies, 
obtained from bare fermion and bare bosonic propagators 
{G°W° scheme). 



principle be well described by the lowest order pertur- 
bative Majorana self-energies (we will see in Sec. 



me 



that this naive statement is not completely correct in 
some specific frequency range). These self-energies are 
obtained at order a by a single boson exchange, see Fig.[8j 
where both fermionic and bosonic lines are given by the 
bare propagators Eqs. ( 22p4 ) and Eq. (31 1, respectively. 
In terms of resummation of diagrams, this level of ap- 
proximation corresponds to the so-called G°W° scheme 
for the generic fermion-boson models arising in the fields 
of heavy fermions^ or strongly interacting electron liq- 
uids^ (G° denotes the bare fermionic propagator and 
W° = G% the bare bosonic one). It is alternatively called 
the (non-self-consistent) Born approximation in the con- 
text of disordered systems,^ and it roughly corresponds 
to the weak-coupling Bloch-Redfield approximation for 
the real-time non-equilibrium dynamics of the spin-boson 
model.^ Note that the present Majorana diagrammatics 
is formulated for the equilibrium dynamics only. Exten- 
sion to temporal dynamics under a sudden switching (of 
the bath or transverse field) could be incorporated in our 
formalism by using two-times Keldysh Majorana Green's 
functions, which lies however beyond the scope of the 
present work. Standard evaluation of the diagrams pro- 



vides the following imaginary part of the Majorana self- 
energies at zero temperature: 



G"w< 



(lu) = a 



-XmG°Je)lmG%(uj - e)(34) 



G"w" 



J de 



(w) = a / ^-XmGl(e)lmG%(u> - e)(35) 



Using the bare propagators ( 22p4 ) and ( 30 1 , one readily 
obtains: 



-7raw c 0[wc " (M - A) 2 }6(u 2 - A ) 
: K (36) 

(37) 



-nau> c 8(uj 2 — to 2 ) 



M 



We are now equipped to start understanding the NRG 
results discussed in Section II D| Owing to the general 
relation Eq. ( 20 1 between spin autocorrelation functions 



and Majorana propagators, we restrict the discussion to 
the latter in the rest of the paper. For a <C 1 one can 
safely neglect the small real part of the self-energy (cor- 
responding to a tiny Knight shift of the spin precession 



frequency A), so that we get from Dyson's equations ( 26 
28| for \uj\ < uj c : 



2mG r . 



lmGf w ° 



ttoiuj 



l-s 



\ui\ - A| s #(u; 2 - A 2 ) (38) 
(39) 



(cj 2 - A 2 ) 2 +w 2 [Tzaul- S \u\ 



lmGf w °(u;) 



A 2 



A 2 U 



(lo 2 - A 2 ) 2 +w 2 [iraujl- s \uY 



(40) 



Following the discussion given in Sec. II D we start 
by interpreting the transverse susceptibilities XmGy w 

and lmG^ 0w0 . Clearly the S(u> — A) peak is replaced 
by a sharp resonance, due to the small but finite life- 
time given by TmYi y (oj ~ A) in both expressions ( 39 ■ 
40 1 . The low-energy behavior is however different in 
those two quantities, as we find XmG^ w 

and ZmGf° w ° oc |w| s , in agreement with the NRG re- 
sults, see the two topmost panels in Figs. |3][5] and the 
related discussion in Sec. |II D| A higher density of low- 
energy modes is indeed obtained in C z than in C y , due to 
the coupling of the bosonic continuum to the spin along 
the z-axis. This leads therefore to the following decay 
of the z-component of the spin autocorrelation function 
C z (t) oc l/t 1+s in the long-time limit. In the perturbative 
regime a<l and a <ti cv c , it is not surprising to find very 
accurate agreement between the NRG and our analytical 
expressions (note that the real parts of the self-energies 
have been included in all plots). Turning to the longi- 
tudinal susceptibility TmG^ w , we observe in Eq. ( 38 ) 
a broad shoulder for \uj\ > A, which is also seen in the 
numerics, but we obtain instead a hard gap in the range 
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FIG. 9: Selfconsistent fermionic self-energies (analogous to 
the GW° approximation), allowing to improve the calcula- 
tion in frequency ranges where lowest order perturbation in- 
correctly predicts the absence of spin excitations. This scheme 
recovers multiparticle excitations that fill in the spurious gaps 
of bare perturbation theory. 



< A, in disagreement with the NRG data. This sur- 



G U W" 



at 



prising failure of the perturbative method for G 
weak dissipation comes from multiparticle effects beyond 
leading order, which do play a qualitative role whenever 
the spectral density is identically zero at lowest order in 
perturbation theory (otherwise higher order self-energy 
corrections are always small at weak coupling). We now 
examine this question in greater detail. 



C. Multiparticle effects 

We still consider here the spin dynamics at very weak 
dissipation a«l and a <C a c , with the goal to improve 
the diagrammatics for the the lowest order longitudinal 
susceptibility ImG^ w , whose spectrum is spuriously 
gapped at lowest order in perturbation theory. This 



inconsistency can be easily tracked to expression (22 1, 
where the bare transverse propagator G° resumes to a 
S(\ui\ — A) peak and misses the uj 2+s low-energy modes 
obtained from the coupling to the bosonic bath, see equa- 
tion (39 1. The discrepancy is thus resolved by reinjecting 
the full Majorana-Green's functions Gi into the Majo- 
rana self-energies (this is the so-called GW° scheme, see 
Fig# 

TmY% w \u) = a f —lmG y (e)lmG%(u) - e)(41) 



de. 



a I —XmGJe)XmG%(u)-e)U2) 
/o t 



The low-energy tails of the transverse susceptibility will 
therefore fill in the gap of the longitudinal susceptibility, 
so that in principle a single iteration of the above self- 
consistent equation should be sufficient required. The 
results are shown in Fig. |10| demonstrating our correct 
physical interpretation of the missing low energy tails 
in G x . We note that the two transverse susceptibilities 
are very mildly affected by the self-consistency at weak 
coupling, and were therefore not shown. 
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FIG. 10: Longitudinal spin susceptibility C x (uj) for A/uj c = 
0.1, s = 0.1 and a = 0.001,0.004,0.007 (bottom to top), 
comparing the NRG result to the self-consistent GW° scheme, 
allowing to recover the low-energy tail below the threshold 
ui — A (results are quantitatively correct only for the smallest 
a value, see Sec. |IVE] for the improvements that are required 
near the quantum critical point). 



A second instance where multiparticle effects are rele- 
vant occurs in the frequency range above the high energy 
cutoff, now for all three spin correlations. Indeed, lowest 
order perturbation theory provides again a gap for ex- 
citations with |w| > u C) see Eqs. ( 36p7 ) in the case of 
one of the transverse spin correlation function (similar re- 
sults are obtained for the other spin susceptibilities). In 
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FIG. 11: Transverse spin susceptibility C z {u)) for A/cj c = 0.1, 
s = 0.1 and a = 0.002, comparing the raw NRG data to the 
self-consistent GW° scheme. In contrast to bare perturbation 
theory, which leads to a gap for cj > ui c , the self-consistent 
resummation correcly describes the tail arising from multi- 
particle excitations above the threshold ui c (note the linear 
frequency scale). 

contrast, the NRG data show a tail with small yet non- 
zero spectral weight above the cutoff ui c , see Fig.[Tl] This 
continuum of magnetic excitations decays very quickly at 
increasing energy, and displays furthermore well-defined 
thresholds, pointing again to multiparticle effects, where 
more and more bosons are exchanged with the spin de- 
grees of freedom, in agreement with recent observation!^ 
Using the fermionic self-consistent scheme proposed in 
Eqs. (41 42), we indeed reproduce numerically both the 



rapidly decaying tail and the multiparticle thresholds at 
quantized values of the cutoff lj c (we note that the multi- 
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particle thresholds seen in the raw NRG data shift do not 
precisely match the quantized values, possibly an artifact 
of the discretization procedure). This final comparison 
concludes our study of the weak coupling regime a -C a c 
of the spin-boson model. 



IV. DYNAMICS NEAR THE QUANTUM 
CRITICAL POINT 



y 



G y 



G° 



G'<i> 



A. Role of the bosonic mass 

We provide here general considerations on the struc- 
ture of perturbation theory in terms of the coupled field 
theory of fermions and bosons, with the emphasis on the 
regime near the quantum critical point. The main moti- 
vation is to understand precisely where the singularities 
associated to critical modes will occur in the diagram- 
matics. From the knowledge gained in the previous weak 
coupling analysis, we understand that the fermionic sec- 
tor is immune to singularities as long as the bosonic prop- 
agator remains regular. Clearly self-consistency for the 
fermionic self-energies (GW° scheme) provides only small 
perturbative correction and cannot account for critical 
behavior. One can convince oneself that the frequency- 
dependent vertex V (see Fig. [7|) is also regular as dissipa- 
tion is increased. The only possibility left is that only the 
bosonic propagator encapsulates the singular dynamics 
via its self-energy (see general diagrammatic expression 
on the lower panel in Fig. [7|. Indeed using the exact 
Dyson's equation ( 29 ) and the asymptotic bare propaga- 
tor ( 32p3 ), we get the low frequency retarded bosonic 
Green's function: 



G*(w) 



1 



m$ + [a + ibSign(uj)]\uj\ s — S$(w) + X$(0) 

(43) 



4w c 



-S*(0) 



(44) 



with m$ the bosonic mass, a and b real coefficients that 
can be read from Eqs. ( 32|[33 |. This simple expression 
allows to understand the development of critical fluc- 
tuations in the spin-boson model: the quantum critical 
point is just associated to a vanishing of the mass m$ in 
the denominator of the bosonic propagator, due to the 
progressive building up of the bosonic self-energy. At 
the quantum critical point, S$(0) = s/(4w c ), and the 
correlation function now diverges as G$(w) oc |w|~ s at 
low frequency, owing to expression ( 43 ) . One can check 



that the frequency dependence of the bosonic self-energy 
E$(w) is subdominant with respect to uj s in all orders in 
perturbation theory, so that the — s exponent of the crit- 
ical fluctuations turns out to be exact. In contrast, the 
bosonic Green's function has a non-divergent power law 
behavior at low frequency G$(o>) — G$(0) oc \uj\ s in the 
whole delocalized phase. This change of exponent near 
the quantum critical point can be used to explain the 
behavior of the transverse spin susceptibility, calculated 



FIG. 12: Fermionic self-energies (without self-consistency) in- 
cluding a dressed bosonic propagator (this is analogous to the 
G°W approximation). This scheme is necessary to recover the 
critical modes near the quantum critical point, but asks for a 



careful evaluation of the bosonic self-energy (see Fig. 13 and 
the related discussion in Sec. IV B and Sec. IV D I. 



previously with the NRG. When using the full bosonic 
propagator within the fermionic self-energy (see Fig. 12 1, 
we get: 



-Sign(w)ImG$(oj) 



cx |cj| s at a = a c 



(45) 
(46) 



From the exact Dyson equation for the r\ z Majorana 
Green's function Eq. ( 28 ) , we then get at low frequency: 



TmG 



G U W 



oc 



y 

at a 



Oir 



(47) 
(48) 



The divergence of the transverse spin susceptibility at the 
quantum critical point is thus explained (again, we stress 
that the above exponent is exact). Considering also the 
second transverse susceptibility Eq. p7| , we get at low 
frequency: 



XmGf w M 



oc \lu\ 



2-s 



at a — a c 



(49) 
(50) 



as also observed with the NRG calculations of Sec. Ill Dl 
This discussion shows that even in the weak dissipa- 
tion regime a <C 1, the presence of such critical modes 
for a < a c will invalidate bare perturbation theory. It is 
thus necessary to investigate in more detail the bosonic 
correlations, in order to improve the agreement between 
numerics and diagrammatic calculation. We start by an 
analytic computation of the bosonic self-energy up to 
two-loop order. 



B. Critical line at two-loop order 

The lowest order (one-loop) contribution to the 
bosonic self-energy (first diagram in Fig. 13 ) is readily 
computed: 



^lloop 



2A A 2 - (uj + iO+y 



(51) 
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FIG. 13: Diagrammatic representation for the bosonic self- 
energy computed up to two-loop order, leading to Eq. (1531). 



This level of approximation in the diagrammatics is 
equivalent to the standard RPA for the interacting 
electron-gas From vanishing of the mass in Eq. (44), 
the above result gives the critical dissipation at one-loop 
order 



Hoop 



sA 

2 LOn 



+ 0(s 2 ), 



(52) 



recovering previous results.^ We recognize here the 
multi-scale nature of the perturbative expansion in the 
dissipation strength. While the high-energy part of the 
spectrum is well described as long as a < 1, the low 
energy sector will be immune to quantum critical fluctu- 
ations as long as a <C a c . While a c vanishes for s — > 0, 
showing that the quantum critical point is perturbatively 
accessible, the above result is only accurate for very small 
s. As we will demonstrate now, an anomalous power-law 
also shows up in the dependence of a c with respect to 
the dimensionless parameter A/w c . To establish this re- 
sult, we need to push the weak-coupling calculation of 
the bosonic mass to next to leading order (two-loops), 
see Fig. [13} After straigthforward but lengthy calcula- 
tions, we find at zero frequency: 



^21oop/ n x _ 
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(A + e) 2 
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(53) 



where the last equation applies for A < cj c . We now solve 
the equation s/(4w c ) = E^ loop (0) + S 2 i loop (0) determining 
the quantum critical point, which gives in the small s 
limit: 



v 21o °P(w) - 



£ — 

2^ 



1 



i~ sl0g (^ 



C(s 3 ). (54) 



This precise form confirms that the weak-coupling dia- 
grammatics is controlled as a systematic small-s expan- 
sion. The presence of a logarithmic correction hints for 
the following anomalous power-law behavior: 



21oop 



+ (£)' 



(55) 



We will first examine this analytical result with respect 
to the NRG phase diagram, before turning to a more 
rigorous derivation. 



Phase diagram of the sub-ohmic spin-boson 
model 



According to the above result, the formula (55) for 
the critical dissipation a c remains perturbative either for 
small s (at arbitrary A/w c < 1) or for small A/uj c . In 
the latter limit, s can take relatively large values, but 
should not be too close to 1, otherwise the small co- 
efficient (A/wc) 1 "* becomes of order 1, and perturba- 
tion theory breaks down. In this case, an alternative 
approach, based on a develop ment a t small ft = A/lo c 
but finite a can be put forward! 19 * 21 1 This results in the 
well-known one-loop scaling equations (which were ini- 
tially obtained from the long-range unidimensional Ising 
model) : 



da 

~dl 
dh 

~dl 



(1 — s)a — ah 2 
(1 - a) ft 



(56) 
(57) 



where the parameter I characterizes the change of scale 
during the renormalization procedure. We stress that 
these flow equations are valid for all values of a, as long 
as h <C 1, and are thus complementary to the regime of 
validity of our small a expansion at fixed ft. Combining 
both equations, one obtains: 



d 
Jl 



log(a) 



ft 2 
2 



(l-s)log(ft) 



= 



(58) 



From Eqs. (56) and (57) the renormalization fixed point 

s for s < 1, so that the 



occurs at a 



1 and ft* 



critical line a c (h) is obtained by the condition: 



log(a c )+(l-a c )- 



h 2 - (1 - s) 



-(l-s)log 



ft 



= o 

(59) 

leading to the behavior a c oc ft 1 s at h <C 1 for all s. We 
also recover from this analysis the exact value a c = 1 for 
s = 1 in the h — > limit. A joint comparison of the small 
a expansion and of the Kosterlitz flow to the NRG data 
is provided in Fig.[l4]for the case h = 0.1. While the two 
expansions work where they are expected to, the two- 
loop small a-development is quantitatively accurate on 
the wide range < s < 0.8, while the one-loop small h- 
development is restricted to the small range 0.9 < s < 1. 
Although pushing the small h expansion to higher order 
would enhance the range of validity of the method, it is 
unlikely that this would recover the expected vanishing 
of a c ~ (s/2)h at small s. In contrast, the small a- 
expansion converges very rapidly even for the largest s 
values: the extreme case s = 1 leads from Eq. (55) to the 
successive estimates ev lloop = 1/2 and cv 21oop = 3/4 in the 
A/uj c <C 1 limit, approaching the exact value a c — 1 in a 
steady fashion. As a clear illustration of this statement, 



we finally present in Fig. 15 a more systematic study of 
the phase diagram for several h values, showing that the 
small a expansion remains accurate even at vanishing ft, 
provided that s < 0.8. 
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FIG. 16: Upper panel: ladder series for the bosonic self-energy 
that is needed to obtain the correct phase boundary within 
the diagrammatic calculation. Lower panel: single frequency 
approximation used at the two-rung level (and similarly for 
higher order contributions), see text. 



FIG. 14: Phase diagram in the (s,a) plane for h = A/co c = 
0.1, comparing the NRG data to the Kosterlitz flow obtained 
from the resolution of Eq. ( 59 1 (giving reliable results for 0.9 < 



s < 1), and to the one-loop and two-loop formulas (521 and 
(55 I (the latter is accurate on the wide range < s < 0.8). 
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FIG. 15: Phase diagram in the (s,a) plane for h = A/co c = 
0.00001,0.0001,0.001,0.01,0.1 (bottom to top), comparing 
the NRG simulations versus the two-loop weak dissipation 
result (valid for s <C 1). 



D. Ladder resummation in the bosonic self-energy 

We have just checked that our analytical expression 
for the critical dissipation ( 55 ) is in excellent agreement 



with the NRG data. However, the derivation of this re- 
sult was performed by a strict expansion at two-loops, 
and a re-exponentiation of the logarithmically singular 
terms. In the spirit of a numerical scheme based only 
on Green's functions, such as developed in Sec. |III| the 
correct phase boundary will not be recovered if strict 
fourth order bosonic self-energies are used, and a clever 
resummation scheme should be found, that would be a 
diagrammatic equivalent to the re-exponentiation of the 
analytical result. Insight on this issue can be gained from 
the analytical derivation, because the logarithmic term 
in Eq. (54 1 can be traced back to a single diagram at 



two-loop order, namely the one with a single-rung ladder 
exchange in Fig. 13 (the last two diagrams with Majo- 
rana self-energy corrections are regular, and provide the 
perturbative corrections of order s 2 in Eq. (55)). This 



quantitative phase diagram, will be obtained upon per- 
forming a ladder resummation to all orders, simila r to 
the diffuson mode in disordered systems,^ see Fig. 16 
Performing such calculation, even numerically, turns out 
to be difficult, because this requires the self-consistent de- 
termination of a four-electron vertex (depending of three 
independent external frequencies). A standard approxi- 
mation in the field of disordered systems, 55 but also in 
the context of interacting electron liquids^ (called the 
Ng-Singwi scheme, or GW c g approximation) amounts to 
keeping a single running frequency within the vertex, 
see lower panel of Fig. [16] in the case of a two-rung 
ladder. In our case, this approximation is vindicated 
by the fact that the low frequency part of the bosonic 
Green's function provides the most important contribu- 
tion to the Feynman diagram, so that strict energy con- 
servation can be omitted from one rung to the next. We 
stress again that the inclusion of such vertex corrections 
in the fermionic self-energies are regular, and can be ne- 
glected at weak coupling, contrary to the singular bosonic 
propagator that we treat here. The ladder series includ- 
ing the full three particle vertex can now be explicitely 
carried out for the zero-temperature Matsubara bosonic 
self-energy: 



bladder / 



(iv)~a / ^G° v (uj + iv)G° x (-icj) 



where 
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de 



1 - A°(iw, iv) 
(60) 



A°(kj,iu) ~q / — G%(ie)G°(-iu-ie)G° x (uj + iv + ie). 
J 2n 

(61) 

Evaluation of the zero-frequency self-energy can be car- 
ried out analytically in the small A/lo c limit, leading to: 
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(62) 



The phase boundary can now be obtained by solv- 
ing E^ dder (0) = l/G%{0), leading to the critical valu< 
gadder = ^ 2 ) ( A/u c ) , in agreement with Eq. ((55 1 



hints that the correct bosonic self-energy, and hence a 



Note that the small corrections at order s are missing 
in the above derivation, because fermionic self-energy in- 
sertions were discarded in the present calculation (these 
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were computed explicitely in Sec. IV B). This correct re- 



sult for the critical dissipation confirms that the ladder 
resummation is the good strategy to reproduce the cor- 
rect phase boundary in a fully diagrammatic calculation. 
We finish the paper by completing this study of the phase 
diagram, which will require a careful discussion of vertex 
corrections in the fermionic self-energies as well. 



E. Full Majorana diagrammatics and the 
three-particle vertex 

We are now equipped to collect all the needed ingre- 
dients for a successful implementation of the Majorana 
diagrammatics up to the quantum critical point (as long 
as a <C 1, i.e. for values of s not too close to the ohmic 
regime) . The next step beyond the previous levels of ap- 
proximation (G°W° and GW° schemes) is the inclusion 
of the vertex corrections (at ladder level) in the bosonic 
self-energy Eq. (62), which is required to obtain the cor- 



rect phase boundary and to describe quantitatively the 
critical low-energy fluctuations. We note that the GW 
scheme with fully self-consistent propagators but with- 
out vertex corrections does not recover the correct criti- 
cal point, and will not be further considered. We there- 
fore propose a GWr approximation, which should incor- 
porate the vertex corrections both for the bosonic and 
fermionic self-energies (see Fig. [7]). Indeed, in order to 
emphasize the importance of the vertex corrections in 
the fermionic sector, let us first rewrite the bosonic self- 
energy as E$ der (u>) ~ a R /(2A), introducing a renor- 
malizcd coupling 



= sfaT 



2otU) c 



(63) 



in terms of the static three-particle vertex V (see Fig. [7]) 
computed at the ladder level. This quantity is plotted 




FIG. 17: Renormalization factor a R /a of the dissipation 
strength from Eq. ( 63 I computed for A/lj c — 0.1 as a function 
of a up to the critical point (a c = 0.007,0.105 respectively 
for the two values s = 0.1,0.5). 

in Fig. [17] for two values of s as a function of a < a c , 
which shows that the full vertex T remains finite up to 



the critical point, in agreement with our general expec- 
tation that this quantity remains regular at all orders in 
perturbation theory. Interestingly, the renormalization 
of the coupling constant a can be relatively large, up to 
20% and 75% for the two cases s = 0.1 and s — 0.5 con- 
sidered here. Such a sizeable effect will impact as well 
the fermionic self-energies, so that our consistent GWr 
scheme will finally rely on the following set of equations: 



' / —XmG v (e)lmG<s,(oj - e) (64) 
Jo T 



n" / — XmG x (e)XmG^{uj — e) (65) 



a R 
2A 



(66) 



with a given in Eq. ( 63 1 . The numerical implementa- 



tion of the above equation, with self-consistency imposed 

ited, 
(we 



by Dyson's equations (|26|29|, is readily implemented, 

and 
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and gives the results displayed in Figs, 
note again that fermionic self-consistency is only used in 
order to fill in the spurious gap of C x , but it is not too 
important here). In the case s — 0.1, the outcome of the 
Majorana diagrammatics is at par with the NRG results 
for all three spin susceptibilities, both at low and high 
energies, and including the quantum critical regime. For 
the larger s = 0.5 value, agreement is still quite good, 
and the visible deviations in the quantum critical regime 
(for a = 0.1) arise from slight discrepancies in the critical 
a c between the diagrammatics and the NRG, so that the 
crossover scale is not accurately captured. We emphasize 
that this quantitative description of the spin dynamics re- 
lies of the weak coupling nature of the spin boson model 
for small s values. The key step here was to understand 
that the most important contributions to the diagram- 
matics are the renormalization (beyond leading order) of 
the bosonic mass shift and the vertex corrections. It is 
possible that analytical renormalization group techniques 
(such as developed in Ref.[57]for the quantum phase tran- 
sition in the pseudogap Kondo problem at weak coupling) 
may be applied to the spin-boson Hamiltonian based on 
the present understanding, but this would likely require 
more complex machinery to obtain similar results. This 
successful comparaison seals our study of the spin dy- 
namics in delocalized phase of the spin-boson model. 



V. CONCLUSION 

We have investigated the equilibrium spin dynamics of 
the sub-ohmic spin-boson model, which generically de- 
scribes a two-level quantum system coupled to a dissipa- 
tive environment of harmonic oscillators (characterized 
by a power-law spectrum J(uj) = 2irauj]r s (jj s 6{<jj)6(uj c — 
w) with s < 1) and submitted to a perpendicular mag- 
netic field A. Two complementary techniques were em- 
ployed, that were quantitatively benchmarked against 
each other. First, the bosonic Numerical Renormaliza- 
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FIG. 18: Spin susceptibilities C z {uj), C y (u)) and C x (uj) (top to 
bottom panels) for A/uj c = 0.1, s — 0.1 and increasing dissi- 
pation strength a = 0.001, 0.004, 0.007 (bottom to top), com- 
pared to the full Majorana diagrammatics (GWT scheme). 



0.001 



FIG. 19: Spin susceptibilities C z {uj), C v (uj) and C' x (lj) (top to 
bottom panels) for A/uj c = 0.1, s — 0.5 and increasing dissi- 
pation strength a = 0.02, 0.06, 0.1 (bottom to top), compared 
to the full Majorana diagrammatics (GWf scheme). 



tion Group (NRG) was used with an optimized broad- 
ening method ("b-trick") to obtain reliably the dynam- 
ical spin susceptibilities from weak (a < 1) to strong 
(a ~ 1) dissipation. Here, we considered not only the 
usual transverse spin susceptibility (orthogonal to the 
field and parallel to the bosonic bath), but also for the 
first time the longitudinal (parallel to the field) and the 
second transverse (orthogonal to both bath and field) 
ones. Second, we established a general Majorana-fermion 
diagrammatic perturbation theory, controlled in the weak 
dissipation limit a -C 1. Because the spin- localization 
quantum phase transition takes place at a critical value 
a c which becomes perturbatively accessible [a c <C 1) ei- 
ther in the small field limit (A <C w c for all s < 1) or 
in the strong sub-ohmic regime (s -C 1 for all A/lu c ), it 
was possible to extend the diagrammatic method up to 
the quantum critical point. We then examined precisely 
how various flavors of GW-like approximations perform 
against the numerically exact results, in two different 



regimes. For a -C a c , bare perturbation theory (anal- 
ogous to the standard Bloch-Redfield approximation for 
the non-equilibrium dynamics) was sufficient to describe 
well the finite energy dissipative features near the ap- 
plied magnetic field A, but fermionic self-consistency was 
needed in order to improve both on the low and high 
energy parts of the spectrum. Self-consistency was re- 
quired because some of the bare self-energies can show 
spurious gaps that are in reality filled by multiparticle 
excitations. For a < a c , huge deviations appeared in 
the low-energy limit between lowest order perturbation 
theory and the NRG data, due to the proximity to the 
quantum critical point. We showed that the effective 
interaction mediated by the bosonic field was the cru- 
cial ingredient to correctly describe the physics close to 
criticality, and that it should include vertex corrections 
(at the ladder level) in the bosonic and fermionic self- 
energies for quantitative agreement with the numerics. 
Such diagrammatic contributions are difficult to evaluate 
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for higher dimen sional field theories of coupled fermions 
and bosons^^J and are rarely considered (see however 
Ref . I5"8"j) . A two- loop analysis allowed us also to obtain 
a simple equation for the critical line, that describes the 
NRG phase diagram very accurately in a wide range of 
parameters. We also presented a proof of Shiba's relation 
generalized to the sub-ohmic case and the derivation of 
an effective bosonic action that may support the validity 
of the quantum-to-classical mapping. 

We finally discuss possible outlooks for the method- 
ology and the ideas developed here. First, the addi- 
tion of a magnetic field component eS z parallel to the 
quantization axis may be considered with special inter- 
est^. Indeed, understanding the magnetization process 
is the next step towards a good description of the lo- 
calized phase of the spin-boson model, which was non 
considered in the present work, and which has led to re- 
cent debateP2H2H Second, further extensions of the model 
to multiple bosonic environments present very interest- 
ing physics in the ohmic case,^' which remains to be in- 
vestigated in greater generality. Third, an important 
consequence of our results is that the long time equi- 
librium dynamics is never accurately captured by lowest 
order perturbation theory in the coupling between the 
spin and the environment (unless dissipation is vanish- 
ingly small). This remark should have deep implications 
for the out-of-equilibrium long-time dynamics in dissipa- 
tive models, which ought to be re-examined beyond the 
Bloch-Redfield regime. Extension of the Majorana dia- 
grammatic method to two-times Keldysh Green's func- 
tion could be also considered in the context of sudden 
quenches in the spin boson model near the quantum crit- 
ical point. Lastly, considering the spin-boson Hamilto- 
nian as a toy model for quantum criticality, the quanti- 
tative understanding achieved here may be relevant for 
field theories of electrons coupled to bosons in higher di- 



mensions. 
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Appendix A: Shiba's relations for the sub-ohmic 
spin-boson model 

The goal of this appendix is to provide a simple deriva- 
tion of Shiba's relation Eq. ( 16 ) generalized to the sub- 
ohmic case. The first step of the proof is to note that the 



linear coupling between the z component of the spin and 
the bosonic continuum, see Hamiltonian Eq. |l]), implies 
an exact relation between the longitudinal spin suscepti- 
bility and the full retarded bosonic Green's function: 

G»(w) = G%(uj) + a[G%{u:)] 2 x z (uj) (Al) 

where the bare bosonic Green's function G%{oS) was in- 
troduced in the previous equation ( 30 1 . This connection 



can be easily established from the use of equations of 
motions, or equivalently by using standard manipulation 
of a path-integral representation of the problem. The 
second step lies in a low-energy analysis of the bosonic 
propagator G$(uj) from the diagrammatic perspective. 
The key point is that the bosonic self-energy E$(w), 
which was defined in Eq. ( 29 1 and expressed diagram- 
matically in Fig. [7j has a regular low-frequency behav- 
ior, owing to the role of the transverse magnetic field 
A in the bare Majorana Green's function G°(w), see 
Eq. (23). As a consequence, all bosonic self-energy di- 



agrams are cut in the long-time limit due to this energy 
gap, so that E$(w) = E$(0) + 0{uj) at low energy. Al- 
though the gap will be in reality filled by low energy 
excitations if one considers the full Majorana propagator 
G y {uj), we have seen in Sec. Ill B that a quickly vanishing 
tail of excitations is obtained, so that the above argument 
persists non-perturbatively (provided that perturbation 
theory is convergent). As a concrete illustration, equa- 
tion (51 1 computed at one- loop order demonstrates the 



regular character of the bosonic self-energy in the small 
frequency limit. Accounting for the static part of the 



self-energy within the explicit Dyson equation (43), and 



neglecting the term E$(cj) — E$(0) at low energy, one 
obtains: 
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(A2) 



where the full mass m$ was given in Eq. (|44j) . Similarly, 
the bare correlation function obeys a similar behavior: 



1 



ib 



irr. 



(A3) 



now involving the bare mass m\ . We finally use the exact 
relation ( Al ) in the low energy limit, which imposes two 



conditions by taking the real and imaginary parts: 
11a 
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Solving this system to eliminate the mass m$ leads to 
the simple expression: 
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asymptotically exact in the small w limit. Using the 
known value for the coefficient —bj (m%) 2 = 2ttu>1~ s from 
Eq. (32), and relation ( fl4|), w e finally obtain the gener- 
alized Shiba relation Eq. (|16[). 



Appendix B: Quantum-to-classical mapping 



We can reformulate the usual quantum-to-classical 



mappings in terms of an effective quantum action in 
imaginary time calculated perturbatively in the coupling 
constant a. Indeed, in the e — > limit, we first notice 
that the linear form of the coupling in the Hamiltonian 
Eq. leads to the exact relation: 



<S*> = 



-<*> 



(Bl) 



which can be derived by the same strategy as discussed in 
AppendixfA) This result implies that the critical behavior 
of the impurity spin magnetization (S" 2 ) can be extracted 
from the knowledge of the bosonic static average. The 
static critical properties can thus be alternatively under- 
stood from the point of view of an effective action for the 
boson $ only, which is obtained by integrating out the 
spin degrees of freedom perturbatively in a. This can 
be done diagrammatically using the Majorana fermions, 
leading after some standard manipulations to the follow- 
ing expression in imaginary frequency: 

S cS m = l J ^( m<s> + AWn\^)\ 2 + J rfr M [$(r)] 4 

(B2) 

where A oc cu c 1 s . The first important parameter here is 
the bosonic mass m$ = s/(4w c ) — a/ (2 A) + (D(a 2 ), which 
agrees with the one-loop result of Eq. (p3lj) , and controls 



the distance to the quantum critical point at leading or- 
der. The second crucial term is the <I> 4 interaction with 
coefficient u = a 2 /(16A 3 ) + 0{a 3 ). 

Within the delocalized phase, i.e. as long as the renor- 
malized mass toJ does not vanish, the bosonic spec- 
tral function obeys the low-frequency behavior G$ (iv) oc 
l/(mf + A\u\ s ) oc 1/mf - A\v\ s /(m$ ) 2 , i.e. the slow 
decay in imaginar y ti me G^{t) oc 1/t 1+s . Owing to 
the exact relation ( Al ), we recover the correct small fre- 



quency behavior of the longitudinal spin susceptibility, 
see Eq. (40). The low-energy sector of the spin-boson 



mode l is thus well captured by the effective bosonic ac- 
tion ( |B2[ ), which turns out to be equivalent to an Ising 
model in imaginary time with long-range spin exchange 
decaying as l/(r — t') 1+s , ex actly as expected from the 
quantum/classical equivalenc e 1 1 8 1 1 9 1. 

Let us finally discuss the low energy behavior at the 
quantum critical point. Scaling analysis relative to the 
massless theory = 0) leads to the bare scaling di- 
mension [u] = 2s — 1. Thus, u is relevant when 1/2 < s 
and irrelevant otherwise. Therefore, in the weakly sub- 
ohmic regime 1/2 < s < 1, the phase transition is de- 
scribed by an interacting fixed poinPSl and the magneti- 
zation (<5 Z ) obeys non-trivial classical exponents, as veri- 
fied in the numerical study of Ref. 1191 On the contrary, in 
the strongly sub-ohmic range < s < 1/2, u is irrelevant 
and mean-field theory should apply, a resu lt th at was 
recently debatedP^HH The effective action (B2) would 



therefore tend to strengthen the idea that the quantum- 
to-classical mapping is robust, unless the Landau param- 
eter u picks up non-analyticities (such as discussed in 
Ref. 1601611 for field theories involving coupled fermions 
and bosons in higher dimensions). However, the regular- 
ity of the diagrammatics would lead us to believe that 
this possibility is unlikely for the spin-boson model. 
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